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Why is it Difficult to Detect a Millisecond Pulsar in Neutron Star 

X-ray Binaries? 

Lev Titarchuk Wei Cui^ and Kent Wood ^ 
ABSTRACT 

We explain why it is possible to detect directly X-ray emission from near the 
surface of the neutron star (NS) in SAX J1808. 4-3658 but not in most other 
low-mass X-ray binaries (LMXBs), with the exception that emission from the 
surface can be seen during bursts events. We show that the X-ray emission from 
SAX J1808. 4-3658 mostly originates in the Comptonization process in a relatively 
optical thin hot region (with an optical depth tq around 4 and temperature is 
around 20 keV). Such a transparent region does not prevent us from detecting 
coherent X-ray pulsation due to hot spots on the NS surface. We give a precise 
model for the loss of modulation, such suppression of the QPO amplitude due 
to scattering can explain the disappearance of kHz QPOs with increasing QPO 
frequency. We also formulate general conditions under which the millisecond 
X-ray pulsation can be detected in LMXBs. We demonstrate that the observed 
soft phase lag of the pulsed emission is a result of the downscattering of the 
hard X-ray photons in the relatively cold material near the NS surface. In the 
framework of this downscattering model we propose a method to determine the 
atmosphere density in that region from soft-lag measurements. 
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1. Introduction 



Low-mass X-ray binaries (LMXB) presumably contain a weakly magnetized neutron 
star (NS). Recently Titarchuk, Bradshaw & Wood (2001), hereafter TBW suggested a new 
method of estimating B- field strength using the magnetoacoustic oscillation model and found 
that the B-field strength for a number of NSs in LMXBs is ~ 10^ gauss. Near-coherent mil- 
lisecond X-ray pulsations have been observed in 4U 1728-34 during thermonuclear (type I) 
X-ray bursts (e.g. Strohmayer et al. 1996). They are interpreted as X-ray intensity modu- 
lated at period close to the spin period of the neutron star. Thus these B-field estimates in 
addition to the detection of the milisecond pulsation give strong arguments for LMXB neu- 
tron stars to be progenitors of millisecond radio pulsars (MLP) (see review by Bhattacharya 
& van den Heuvel 1991). But there is still a question why these coherent pulsations are not 
found in persistent emission despite careful searches (Wood et al. 1991; Vaughan et al. 1994). 
The lack of coherent pulsations has been explained as modulation loss from gravitational 
lensing (Wood, Ftaclas & Kearney 1988; Meszaros, Riffert & Berthiaume 1988) or from to 
scattering (e.g. Bainerd & Lamb 1987; Kylafis & Khmmis 1987 ). The third explanation 
for the lack of pulsations is presented by Gumming, Zweibel & Bildsten (2001) who argue 
that the surface field is weak because of magnetic screening ^. In this Letter we put forth 
arguments for smearing out of the pulsar signal due to electron scattering in the optically 
thick environment typical of most observed LMXBs. In §2 we study the scattering effect and 
its relation with the observed timing and spectral characteristics for various QPO sources. 
In §3 we investigate the scattering effects in QPO sources. In §4 we present results of the 
RXTE data analysis of spectral properties of X-ray radiation in SAX 1808.4-3656. In §5 we 
analyze the downscattering model and its application to the observed soft lag phenomenon 
detected in the coherent pulse signal from SAX 1808.4-3656. Conclusions follow in §6. 



2. Scattering Effect on QPO Pulse Profile and Pulsar Profile 

The transfer problem to be solved is that of a sinusoid, representing either (a) an 
intrinsic pulsar light curve or (b) a high frequency QPO, but in each case seen through a 
surrounding medium that reduces the pulse modulation by scattering. The signal whose 
modulation is affected lies in the range of a few hundred Hz to perhaps 1 kHz. We thus must 
analyze smearing out of a signal at very high frequencies and investigate the dependence 



"The recent discoveries of MLP in XTE J1751-305 (Markwardt ct al. 2002) and in XTE J0929-314 
(Galloway et al. 2002 ) with extremely low mass transfer rates support the suggestion of the absence of the 
magnetic screening for the low mass accretion rates. 
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upon both the intrinsic frequency and the optical depth of the scattering environment. The 

geometry used derives from earher works (see e.g. Sunyaev & Titarchuk 1980, hereafter 
ST80; Titarchuk, Lapidus & Mushmov 1998, hereafter TLM98) modehng energy spectra 
and power density spectra (PDS) of LMXB sources. For the high frequency QPO, [case (b)] 
the TLM98 model identifies the site where the QPO is generated^ which will be at radius 
substantially larger than that of the star. Thus the effective optical depth of interest in case 
(b) will not be the full optical depth to the star surface but rather the reduced optical depth 
appropriate to that site location. 

One more consideration completes the problem specification. It is well established from 
observation that the kHz frequency QPOs in LMXBs are not fixed in frequency but increase 
or decrease in manner correlated with changes in luminosity, which in turn is believed to vary 
directly as mass accretion rate, hence as source brightens the OPO site moves both inward in 
radius and upward in frequency. Consequently, the total change in QPO power as the source 
varies will be a combination of two effects, namely (i) the intrinsic change in QPO power 
that results simply from physics of shifting the resonance frequency that gives the QPO 
and (ii) the further reduction in power that results from its now being observed at higher 
frequency and through greater optical depths both of which enhance the effects of scattering. 
To model the entire process we therefore require a model for the intrinsic change in power. 
This is provided in the cited transition layer (TL) model (TLM98) by approximating the 
QPO process itself as a damped harmonic oscillator under the influence of a driving force. 
Then, the interaction leads to resonance oscillations for which the Q-value {^qpo/ ^^qpo) 
can be high. The resonance amphtude varies as the product of the amphtude of the driving 
force Aciriv) and 1/vqpo (see Landau & Lifshitz 1965). The default intrinsic model, in the 
absence of scattering, is thus that the QPO amplitude falls off as A{v) = Adriv) / vqpo^ but 
scattering contributes additional reduction. In a source such as GX 340+0, where a pair of 
high frequency QPOs is seen, it is possible to test out these ideas in detail, using the model 
derived below and applying it to both QPOs, keeping the optical depth tq the same for 
both kHz QPOs because they are co-located. Detection of the NS spin period is the crucial 
further issue at stake here. The typical problem is the following: since the NS spin gives 
a coherent periodic signal at nearly constant frequency it should be at relative advantage 
in detectability in comparison with QPOs, yet is not see while the QPOs are seen. The 
key, obviously, is that the QPO sees the reduced optical depth while the modulation at the 
spin period experiences the unreduced optical depth. The occasional exceptions where the 
stellar spin frequency does prove detectable (such as MLPs, plus the transient detections in 
X-ray bursters) should be explicable in terms of circumstances that reduce scattering effects 
in those particular cases below values that generally obtain; moreover those exceptional 
reductions should be confirmed by spectra indicative of the lower optical depth and other 
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supporting context observations. 

For a quantitative model of this problem scattering is handled as a diffusion problem in 
a Green's function [G{t)] treatment. The input pulse is x{t) = Asm{ujt + ipo) where u; — 27ru 
is the NS spin rotational frequency cu or the QPO centroid frequency depending on whether 
it case (a) or (b) and (po is an initial phase. The resulting pulse affected by the scattering is 
z{t) and is given by a convolution of x{t) and G{t): 

z{t) = [ x{t')G{t-t')dt'. (1) 

For times t much greater than the scattering time in the Compton cloud Rt/c the initial 
pulse emission x{t) = A{u)) exp [{ujt + ifoji] is calculated as follows 

poo 

z{t)^ A{uj)eyLY>[i{ujt + ipo)\I{uj), and i"H = / e-^''G{u)du, (2) 

Jo 

ST80 analyzed the response of the scattering medium (the Green's function) for various 
source distributions. They show that for a QPO source embedded in the center of the 
cloud (r = 0) the response of the scattering medium (the Green's function), G{t) can be 
approximated by 

G{t) oc {Ci + Cat"^/^) exp[-3i?ro/4ct - 7r'^CTot/3R{To + 2/3)^], (3) 

where tq, R are optical depth and radius of the cloud respectively, c is the light speed. 
Coefficients Ci and C2 can be estimated as follows (see also Wood et al. 2001) Ci = 
57rV3(ro + 2/3)^, and C2 = T^{'i/Tif'^Tl'^{R/cf^ /A. Using formula (3) the integral I{u) 
is analytically calculated by the steepest descent method (see also Prudnikov, Bruchkov 
& Marichev 1981, formula 2.3.16). We can show that an amplitude B{ijj) of z{t) can be 
obtained as 

B{u) = A{u)T{X) exp {-2X) (4) 

where X{uj,to) = y/[7rTo/2{To + 2/3)]^ + {3Rtou/4c)^ and = d^X^/'^ + d^X-^'^ + do 

is a weak function of X which is a constant and the product of 1F{X) exp {—2X) is almost 
unity for uj much less than the inverse of crossing time ^ = (rro/c)"^. The amplitude B{uj) 
decays exponentially when uj increases. 

3. Consequencies of scattering effects for QPO observations: Evidence of the 

resonance excitation in the QPO sources 

As it is seen from formula (4) the amplitude of the NS pulsations with a; ~ 27r x 400 
Hz decreases very rapidly with tq for 

To ^ 4[(c/i?)/5 X 10^ s-^]/[a;/(27r X 400 Hz)]. (5) 
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A radius R of 60 km was chosen as a typical radius of the Compton cloud derived from 
the observations in the frameworks of the TL model (see TBW). In Figure 2 we present the 
rms amplitude of the high frequency QPOs as a function of kHz frequencies for the Z-source 
GX 340+0 using the results of Jonkcr ct al. (2000). We compare them with the theoretical 
model presented by formula (4). The amplitude of the lower kHz QPOs is perfectly fitted 
by l/uj dependence (see the dashed line on the left side of Fig. 2). Xred — 0-^3 for this 
fit, but the deviation of B{u;) due to scattering (the solid fine on the right side of Fig. 
2) from the 1/u; law (the dashed line on the right side of Fig. 2) is clearly seen for the 
higher kHz QPOs. The best fit to the data {Xred — 1-1) by formula (5) with A{u;) oc 1/u; 
is obtained for tq = 2.7 with R/c — 10~^ s. For tq — 2.7 the rms, B{u;) calculated for low 
kHz frequencies is almost identical to l/oo dependence {xled ~ 0.85). Thus we found that 
the driving force amplitude A^r is practically independent of frequency for kHz QPOs. One 
can argue that the driving force power for kHz QPOs can be determined by the observed 
PDS at kHz frequencies. The PDS power law index (3 typically varies around 1 (see details 
in Wijnands & van der Klis 1998). The resulting resonance amplitude would be uj~^^+l^/^^ 
with an assumption that this power continuum is a driving force reservoir for the resonance. 
As an example we fit this model with /3 = 1 to the data for the lower kHz QPOs 

only, where the scattering effects can be neglected (see inequality 5). We found that for this 
model Xred ~ 1-37. The rms-frequency dependence for the lower kHz QPOs is not affected 
by scattering; it affects mainly the higher kHz QPOs (see Eq. 5). The scattering effects 
make the rms-frequency dependence much steeper for the higher kHz QPOs which does not 
allow us to extract the true driving force amplitude from the data (clearly, it is not feasible 
to dispense with scattering and fit both QPOs solely with type of driving force model) ^. 
It is well estabhshed in the QPO observations that there is a correlation between the QPO 
frequency and the mass accretion (or the count rate). The high QPO frequency increases 
with the count rate reaching a certain level and then it disappears as the low frequencies 
are still detected (Zhang et al. 1998; Mendez et al. 1999). This effect can be explained by 
the exponential decay of the QPO amplitude with increase of tqcj which is a result of the 
simultanenous rise of tq and u with the mass accretion rate, (see Eq. 4). The high frequency 
QPOs presumably originate in the outskirts of the Compton cloud and the optical depth of 
the QPO source location is relatively small (1-3) even in systems with high accretion rate. In 
contrast, the source of the NS pulsations is deeply embedded in the Compton cloud, whose 
For high mass accretion it is hardly to expect the detection of the pulsations from the NS 
surface. The Compton cloud optical depth can be very large. The NS pulsations would be 



^Fits of the rms amplitude vs frequency similar to that presented in Fig. 2 for GX 340+0, have been 
produced for GX 5-1, Cyg X-2, GX 17-2 and Sco X-1. Thus one can conclude that in all these kHz QPO 
observations the rms resonance amplitude 1/ui vs frequency is affected by the photon scattering. 
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totally smeared out by the scattering if the pulsation frequency is about 300-400 Hz and 
higher. In the next section we show that in SAX J1808. 4-3658 the optical depth of spherical 
Compton cloud is relatively low allowing detections. 

4. X-ray spectra of SAX J1808.4-3658 and of Other LMXBs 

We revisit some RXTE spectral data for SAX J1808.4-3658 (Gilfanov et al. 1998; 
Heindl & Smith 1998). The spectrum (which was constructed using data from 30411-01-06- 
00) is well represented by the Comptonization model (CompTT in XSPEC, with spherical 
geometry; Titarchuk 1994) along with some contribution from a blackbody component and 
a Fe Ka— line. The best-fit parameter of the model arc tq = A.5^qIj, electron temperature 
kTe = 19.5l3 2 kcV, Ka energy is at 6.5 ± 0.1 keV. The temperature of the seed photons for 
the Comptonization is within /cTq = O.IIq^ keV. The temperature of the blackbody (BB) 
radiation (presumably from the neutron star surface) is kTj,!, — 0.69 ± 0.02 keV. The quality 
of the fit is high with — 294/322 = 0.91 ^. It is worth noting the relatively small optical 
depth T = 4.5. This can be the main reason why the NS pulsation is detected from this source 
(see inequality 5). If this tq is compared with that for other LMXBs sources (showing QPO 
features) one can find that optical depths for a number of sources arc higher than that for 
SAX J1808.4-3658. They arc 10-11 for Cyg X-2 (Kuznetsov 2002), Sco X-1 and GX 340+0, Tq 
is around 6 for 4U 1728-34 (TBW). On average the seed photon temperature value Tq <^ T^f,. 
Presumably the photon supply for Gomptonization mostly comes from the relatively cold disk 
and the blackbody emission comes from the NS surface. The Gomptonization spectra of the 
QPO sources where the NS pulsations are not observed is indirect evidence for this, van der 
Klis (2000) argues that there is a striking similarity between the power spectra of black holes, 
atoll sources, Z sources and the MLP SAX J1808. 4-3658. This would exclude any spectral 
formation models requiring a material surface and would essentially imply the phenomena 
are generated in the accretion disk around any low-magnetic field compact object. But in 
the case of SAX J1808. 4-3658 we definitely see coherent pulsations of the hard radiation 
which originates in the NS surface [Gui, Morgan & Titarchuk (1998), hereafter GMT98] and 
also we should see the downscattering effect of the hard Gomptonized photons when they 



'^Gierlifiski, Done & Barret (2002), hereafter GDB, also analyzed the X-ray spectrum of SAX J1808.4-3658 
and they showed the spectrum can be fitted by a three-component model (BB component, Comptonization 
component and the Compton reflection). In general, spectral modeling is not unique and it is not by chance 
that our spectral model and global picture are different from that of the GDB. Because at kT^h around 0.7 
keV the iron is highly ionized we conclude that the observed Ka (see §4 and GDB) is formed as a result of 
the disk reflection rather than as the NS reflection. 
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hit the relatively cold material of the NS star (see Fig.l). 

5. Downscattering effect 

For SAX J1808. 4-3658, the observed characteristics of the pulsed X-ray emission (pulse 
profiles and phase lags) and the overall energy spectrum provide useful insight into X-ray 
production processes and the emission environment. Important information regarding the 
density of the NS atmosphere can be extracted from the time lags between different energy 
bands. The soft lags have been detected for SAX J1808.4-3658 by CMT98^. The soft phase 
lag can be due to reflection (Compton downscattering) of the pulsed hard radiation from 
the NS atmosphere. In fact, this effect is unavoidable when hot plasma streams converges 
towards the NS magnetic poles. The hard radiation at the pole illuminates adjoining regions 
where a significant fraction of the hard radiation is reflected. For simplicity we assume the 
input hard photons are monochromatic with energy E^. The photons that emerge from 
the NS atmosphere with lower energy Ei arrive at distant observer later than those with 
higher energy Eh- The delay in arrival time, 5t, is given by ul/c, where u is the number of 
scatterings undergone by the seed hard photon with energy Eh before emerging with energy 
El and / = l/ngO-T is the Thomson mean free path for the photon. The average fractional 
energy loss for the photon of energy E after each scattering is nearly independent of T^;, and 
it is given by 5E/ E ^ —{E — AkThh) / mgC^ ^ —E/meC^ (because 4:kTi,i, <^ E), where mgC^ is 
the electron rest mass. Integrating over multiple scatterings, we have u = meC^{l/Ei — 1/Eh) 
(ST80). The measured soft lag scales as E~^ as shown in Figure 3, and levels off above 10 
keV. In Figure 3 we also show the best fit to the data using the downscattering model 

E-^ ^E^^ + D- St. (6) 

We fit the theoretical E"^ as a function of St to the measured quantities. The best-fit 
parameters are Eh = 20 keV and D = 1.3 x 10"=^ fis'^ keV"^ for which Xred = 11-3/9 = 1.2. 
Because D = UeCaT/meC^ we find that the NS density near the surface is rig = DrUeC^ / car = 
3.3 X 10^^ cm~^. The measured soft-lags in the frameworks of the downscattering models allow 
us to estimate the optical depth r up to what the hard photon can penetrate illuminating 
the NS atmosphere. It takes more than 70 scatterings for 20 kcV photon to reach 5 keV 
energy band. As it is known from the diffusion theory (see e.g. ST80) a photon in the 
course of the random diffusion propagation to optical depth r undergoes on average Ngc ~ 
scatterings. Thus using Nsc — 70 (the scattering number of the input photons which 



''Galloway et al. 2002 also found the hard X-ray pulses arrived up to 770 ytts earlier the soft X-ray pulses 
in XTE J0929-314. 
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finally emerge) one can deduce the penetration optical depth Tpr — ■\/W^j2 — 6. In other 
words the thickness of the NS reflection slab is approximately {mp/aT)Tpr — 15 grams cm""^. 
Photoelectric absorption is suppressed at kT^b = 0.7 keV and thus one can estimate Tpr 
neglecting the photoelectic extinction. The deduced Tpr leads to estimate of the reflection 
albedo A = 1 — (3x^^/4 + 1)^^ = 0.82 (see e.g. Sobolev 1975). It means that more than 80% 
of the input hard photons are reflected by the NS atmosphere. 



6. Conclusions 

Our analysis of the timing and spectral properties of milisecond pulsar SAX J1808. 4-3658 
and comparison of them with those for LMXB QPO sources leads us to these conclusions: 
(1) the detection of the NS pulsations from SAX J1808. 4-3658 was possible because of the 
relatively transparent Compton cloud covering in this source (tq ~ 4). For the majority of 
the analyzed LMXBs the Compton cloud optical depth is at least twice as high as that in 
SAX J1808. 4-3658. The high frequency pulsations with frequencies ^300 Hz are strongly 
wiped out by scattering in the clouds with tq > 4. SAX J1808. 4-3658 is a limiting case 
for this detection. There is a possibihty of flnding the NS pulsations in the sources that 
are much less luminous (than these bright QPO LMXBs), because in them the Compton 
cloud is more transparent (to^4) for the NS pulsed radiation. (2) The kHz QPOs rms vs 
frequency follows the resonance law 1/uo are weakened by scattering in the Compton cloud. 
(3) The soft-lag measurements along with the implementation of the downscattering model 
can provide a tool for density determination near the NS surface. 

We appreciate the fruitful discussions with Paul Ray, Sergey Kuznetsov and particularly, 
with the referee. 
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Fig. 1. — Schematic picture of the model geometry SAX J1808. 4-3658. The hot flow forming 
at the outer boundary of the transition layer streams towards neutron star along magnetic 
field lines. The refiection of the hard radiation of the flow from the neutron star surface 
results in soft time lags of the pulsed emission. 
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Fig. 2. — Examples of RMS amplitude vs QPO frequency for lower and higher kHz peaks 
for GX 340+0 (Jonker et al. 2000). Two solid hues (for lower and higher peaks respectively) 
are theoretical curves for which scattering effects are taken into account. Dashed lines are 
for 1/v law and dash-dotted hne is for law. 
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Fig. 3. — Measured soft X-ray lags with respect to 20-30 keV band along with the best-fit 
downscattering model. 



